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Abstract 

We introduce a nonequilibrium percolation model which shows a self- 
organized critical (SOC) state and several periodic states. In the SOC state, 
the correlation length diverges slower than the system size, and the corre- 
sponding exponent depends non universally on the parameter of the model. 

The periodic states contain an infinite cluster covering only part of the system. 
PACS numbers: 05.40.+j, 05.70.Jk, 05.70.Ln 
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During the past years, systems which exhibit self-organized criticality (SOC) have at- 
tracted much attention, since they might explain part of the abundance of l//-noise and 
fractal structures in nature Their common features are slow driving or energy input 
(e.g. dropping of sand grains, increase of strain, growing of trees) and rare dissipation 
events which are instantaneous on the time scale of driving (e.g. avalanches, earthquakes, 
fires). In the stationary state, the size distribution of dissipation events obeys a power law, 
irrespective of initial conditions and without the need to fine-tune parameters (see for 
examples). In general, essentially three different types of avalanche size distributions are 
possible in systems with slow energy input and instantaneous relaxation events: Either re- 
laxation events occur quite frequently but release only little energy (localized avalanches), 
or they occur very seldom and release a finite portion of the systems energy (avalanches 
which cover the whole system), or the system is at the borderline between both regimes 
and exhibits a power-law size distribution of avalanches. In SOC systems, some mechanism 
pushes the system naturally to the critical point, e.g. the local conservation of sand in the 
sandpile model, or separation of time scales in the forest-fire model. A model which has the 
two noncritical regimes, separated by a critical point, has been introduced e.g. in ||. 

In this paper, we present a model that belongs to the mentioned class of systems, but has 
surprising features not known so far: The region of small avalanches and the region of infinite 
avalanches are separated by a finite region of critical behavior, where the correlation length 
diverges slower than the system size. Since criticality is not restricted to a single point, 
this model can be counted among the SOC systems. In the region of infinite avalanches, 
the system shows synchronization with a period that depends on the value of the control 
parameter. 

Our model can be obtained from the SOC forest-fire model 0] by taking the tree density p 
as the control parameter instead of the ratio between lightning probability and tree growth 
probability f/p. It is defined as follows: Initially, the fraction p of all sites in a square 
lattice with L 2 sites are randomly chosen to be occupied. We call an occupied site "tree". 
The remaining sites are empty. Then we repeat the following rules: (i) Lightning strikes 
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an arbitrary site in the system. If the site is occupied, the whole cluster of s trees which 
is connected to this site (by nearest-neighbor coupling) burns down, i.e. the trees of that 
cluster turn to empty sites, (ii) We grow s new trees at randomly chosen empty sites. 

With these rules, the tree density p is a globally conserved quantity. The model therefore 
does not apply to any kind of burning or reaction process. It would be more realistic to 
speak of "explosions" instead of fires: At a randomly chosen site, an explosion takes place, 
and the whole cluster connected to the explosion site is blown into the air and settles down 
somewhere else in the system. Or one might think of colonies of animals which are dispersed 
into all directions by some enemy or other event. Or, from a purely mathematical point 
of view, one has a nonequilibrium percolation problem. But since the model originally was 
derived from a forest-fire model, we continue to describe it in terms of trees and fires. 

In the following, we discuss the properties of the stationary state as function of the tree 
density. Let s be the average number of trees destroyed by a lightning stroke. For small tree 
densities, there exist only small forest clusters, and consequently only few trees are destroyed 
by a lightning stroke, i.e. s is small. With increasing tree density, the size of the largest 
forest cluster increases, and s increases, too. There exists a critical density pi ~ 0.41, where 
the size of the largest cluster and s both diverge. The critical behavior close to pi can be 
described by exponents which are defined as in percolation theory ||. The size distribution 
of forest clusters is 



where R is the radius of gyration of a cluster. The correlation length is given by £ oc 
(pi ~~ P) u ■ The quantity s can be considered as susceptibility and diverges as s oc (pi — p)~ 7 . 
The exponents are related via the scaling relations 1/er = 7/(3 — r) = pv. In our simulations, 
we found r = 2.15(3), p = 1.96(2), v = 1.20(5), and 7 = 2.09(5). 




(1) 



with a cutoff function C. The fractal dimension of the clusters is defined via 




(2) 
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These exponents are identical to those of the SOC forest-fire model for f/p = (1 — p)/sp, 
if they are appropriately redefined (see |6|-|9[]). In a system which is much larger than the 
correlation length, the difference between a globally conserved tree density and a tree density 
which is only conserved on an average cannot be seen on length scales comparable to the 
correlation length. 

In the region p > p\, which was inaccessible in the SOC forest-fire model, one might 
naively expect an infinite forest cluster which spans the whole system, as in percolation 
theory. However, the values of the critical exponents for p % p\ show that this cannot be 
the case here. The hyperscaling relation d = p(r — 1) is violated |6|||, which means that not 
every part of the system contains a spanning cluster |J. This can be seen in Fig. |](a), where 
in addition to large forest clusters also large clusters of empty sites exist. In contrast to 
ordinary percolation, there is no homogeneously distributed set of large clusters that could 
join at p = p\ to form the infinite cluster. 

Nevertheless, the largest forest cluster has to be infinitely large for p > p\. For p\< p < 
p\ ~ 0.435, we observe the following behavior: The size of the largest cluster s max diverges, 
but slower than L 2 , and the correlation length diverges slower than L. We find s max oc Z/ 1 , 
£ oc L^ 2 , and s oc L^ 3 with p-dependent exponents 01,2,3, while r and p remain unchanged. 
Tab. | shows the values of the exponents for different densities. Fig. |l|(b) shows a snapshot 
of the system for p = 0.43, and Fig. ||] shows the size distribution of fires sn(s) for different 
system sizes at fixed p (a) before and (b) after rescaling. With Eqs. (0),© one can derive 
the scaling relations 0i = p<p2 arid 03 = (3 — r)0i, which are confirmed by the simulations. 
For p — > pi, 01,2,3 approach the values p, 1, and (3 — r)p. 

At p = p 2 ., the correlation length becomes proportional to the system size, and the 
critical state becomes unstable. We observe a discontinuous phase transition to a state 
where the largest cluster contains a finite portion of all trees in the system. Fig. |l](c) shows 
a snapshot of the system for p = 0.45. The system has five homogeneous and equally large 
subphases with different densities. The subphase with highest density contains an infinite 
cluster. When lightning strikes one of the small clusters in this state, only few trees burn 
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down, and the state of the system essentially remains unchanged. When lightning strikes 
the infinite cluster, a large portion of all trees in the subphase with the highest density burn 
down and are regrown at empty sites all over the system. The subphase which used to have 
the highest density now has the lowest density, while the density of the other subphases has 
increased. The values of the five densities are the same as before, except that they are now 
associated with different subphases. If we measure time in units of large fires, the state of 
the system is periodic with a period 5. Increasing p further, one finds four, three (Fig. |l](d)), 
two (Fig. |l|(e)), and finally one (sub)phase (Fig. |l|(f)), where the infinite cluster spans the 
whole system. The shape of the subphases depends on the boundary conditions. 

To understand the occurrence of subphases with different densities, we consider first a 
system with a density far above the percolation threshold 0.59 for random site percolation 
||, and we start with a random initial state. In the first iteration step, lightning strikes 
either the infinite cluster, which consists of nearly all trees in the system, or it strikes one of 
the finite clusters, which are very small. In the first case, only a few small clusters survive 
the fire, and most of the trees are redistributed randomly in the system. In the second case, 
only a small number of trees are redistributed. In both cases, however, the state of the 
system changes little and the new state is close to a completely random state. This remains 
true even after many iterations, since the clusters which survived the first large fire will burn 
down during one of the following large fires, so that correlations cannot increase with time. 

If we decrease p, the clusters which survive the first large fire become larger. When the 
burnt trees are regrown, part of them grow in or near these surviving clusters, where the 
density then will be larger than the mean density. In the space between these clusters the 
density consequently is lower than the mean density. If p is below a threshold p*, this density 
between the surviving clusters becomes smaller than the percolation threshold 0.59. Then 
there exists no infinite cluster in the system after the first time step, and the state with one 
phase becomes unstable. Evidently p* > 0.59. In our simulations, we find p* ~ 0.625. 

For p < p*, a stationary state with two subphases occurs, the density of one subphase 
being above, the density of the other subphase being far below the threshold p*. Decreasing 
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p further, eventually the density of the high-density subphase drops below p*, and a state 
with three subphases is formed, and so on. 

Let pi, . . . ,p n be the densities in a state with n subphases, starting with the highest 
density. Additionally we define the density p n +i of the subphase which contained the infinite 
cluster, immediately after the infinite cluster has been removed from the system, and before 
the removed trees are regrown. As consequence of a large fire, the different subphases just 
exchange their densities, i.e. 

Pi-l = Pi + (Pi - Pn+l) ■ (1 - Pi) I (n (1 - p) + Pl - Pn+l) 

for % = 2, . . . ,n + 1. The last factor on the r.h.s. represents the fraction of trees of the 
infinite cluster that are regrown in the subphase with density pi. We finally obtain 

1 ~ Pi = 1 ~ 92 ^_ = 1 ~ Pn 
1 - p 2 1-/03 " " 1 - Pn+l ' 

Together with p = (1/n) Y%=i pi we have n equations for n + 1 densities. 

In our simulations, we observed a maximum number of n = 5 subphases. In the following, 
we argue that there exists an upper limit for the number of subphases even in an infinitely 
large system: A state with n subphases is only stable if the subphase with the second 
highest density has no infinite cluster. In the limit n —> oo, we would have pi — p<i — > and 
Pi — > p*. In such a state, the density of a given subphase would increase continuously in 
time, until it reaches p*. Then the infinite cluster in this subphase would be destroyed, and 
the whole cycle would restart. We simulated this dynamics, and measured the percentage 
of trees in the largest cluster as function of the density (which itself is a function of time). 
This percentage became finite at p ~ 0.595 < p*, indicating that subphases with densities 
above 0.595 contain an infinite cluster. Consequently states with p 2 > 0.595 cannot exist, 
which gives with Eq. (|3|) a maximum possible number of n — 11 (±2) subphases, and a 
corresponding minimum mean density p = 42(±0.3)%. 

Since the transition between states with different number of subphases is discontinuous, 
there are hysteresis effects. When the density p is decreased, a state becomes unstable 
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when pi falls below p* ~ 0.625, and its number of subphases increases by one. When p 
is increased, however, the subphase with the highest density never becomes unstable, and 
the transition from an (n + l)-subphase state to an n-subphase state takes place when the 
subphase with second highest density P2 starts to contain an infinite cluster, i.e. when p2 
approaches ~ 0.595. 

To conclude, we have described a nonequilibrium percolation model which shows several 
new phenomena which are unknown in equilibrium percolation. Besides the earthquake 
model [|l(]], this is the first model with a SOC state where the correlation length diverges 
slower than the system size and the corresponding exponent depends continuously on the 
parameter. In contrast to the earthquake model, the exponent which characterizes the size 
distribution of avalanches remains constant. In addition to the SOC state, the system shows 
synchronization leading to periodic states with 1 to 5 subphases. 

This work was supported by the Deutsche Forschungsgemeinschaft (DFG) under Con- 
tracts No Dr 300/1-1 and No Schw 348/7-1. 
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FIGURES 

FIG. 1. Stationary states for different densities and boundary conditions (trees are black, empty 
sites are white): (a) p = 0.41, L = 4096, periodic b.c. (b) p = 0.43, L = 4096, p.b.c. (c) p = 0.45, 
L = 4096, p.b.c. (d) p = 0.50, L = 2048, absorbing b.c. (e) p = 0.55, L = 1024, a.b.c. (f) p = 0.63, 
L = 1024, a.b.c. 

FIG. 2. Normalized size distribution of fires for p = 0.43 and L = 512, 1024, 2048, 4096, (a) 
before and (b) after rescaling. 
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TABLES 



p 0.41 0.42 0.43 0.435 

01 1.46(8) 1.72(6) 1.86(8) 1.92(8) 

<p 2 0.79(2) 0.92(3) 0.97(3) 0.99(2) 

03 1.23(3) 1.50(3) 1.62(3) 1.69(5) 

TABLE I. The exponents ^1,2,3 of s max , £, and s for various densities p\ < p < p\ (L = 512, 
1024, 2048, 4096) 
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